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Abstract. An adiabatic, nonradial pulsation study of GW Vir stars is presented. The pulsation calculations are 
based on PG1159 evolutionary sequences with different stellar massess artificially derived from a full evolutionary 
sequence of 0.5895 Mq that has been computed taking into account the evolutionary history of the progenitor 
star. The artificial sequences were constructed by appropriately scaling the stellar mass of the 0. 5895- sequence 
well before the models reach the low-luminosity, high-gravity stage of the GW Vir domain. We compute g-mode 
pulsation periods appropriate to GW Vir variable stars. The implications for the mode-trapping properties of our 
PG 1159 models are discussed at length. We found that the mode-trapping features characterizing our PG 1159 
models are mostly fixed by the stepped shape of the core chemical profile left by prior convective overshooting. This 
is particularly true at least for the range of periods observed in GW Vir stars. In addition, we make asteroseismic 
inferences about the internal structure of the GW Vir stars PG 1159-035, PG 2131+066, PG 17074-427 and PG 
0122-f200. 
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1. Introduction 

GW Vir (or DOV) stars constitute currently one of the 
most exciting class of variable stars, since they represent 
an evolutionary connection between the cool, very lumi- 
nous, asymptotic giant branch (AGB) stars and the hot, 
extremely compact white dwarf (WD) stars. Since the dis- 
covery of their prototype, PG 1159-035, by McGraw et al. 
(1979), GW Vir stars have been the focus of numerous 
observational and theoretical efforts. These stars are be- 
lieved to be the variable members of the PG 1159 family, 
a spectral class of stars characterized by a rather abnor- 
mal surface chemical composition, with atmospheres de- 
void of hydrogen and, instead, rich in helium (~ 42%), 
carbon (~ 43%) and oxygen (~ 15%) (Werner et al. 
1997). PG 1159 stars (80 000 < T^s < 180 000 K and 
5.5 <^ log .9 <^ 8.0) are thought to be the descendants of 
post- AGB stars that, after experiencing a late thermal 
pulse at the beginning of the WD cooling track, return to 
AGB and finally evolve into the hot central stars of plane- 
tary nebulae as hydrogen-deficient pre- WD stars — the so- 
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called "born-again" scenario (Fujimoto 1977; Schonberner 
1979, Iben et al. 1983). In fact, recent evolutionary calcu- 
lations of the born-again scenario incorporating convective 
overshooting have been successful in reproducing the ob- 
served photospheric composition of PG 1159 stars (Herwig 
et al. 1999; Althaus et al. 2005). 

GW Vir variables are low-amplitude, multiperiodic g- 
mode pulsators, with periods in the range from 5 to 30 
minutes. Eleven of such variables are presently known, six 
of which associated with a planetary nebula and some- 
times termed PNNV (Planetary Nebula Nuclei Variable). 
The remainder objects, lacking a surrounding nebulae, 
are called "naked" GW Vir stars. In Fig. ^ we plot the 
observed location of the known GW Vir stars in the 
log Toff — logg plane. In addition, we have included in 
the figure some non- variable PG 1159 stars as well as PG 
1159 stars for which photometric data are not available. 
Note that GW Vir stars are separated into two subgroups, 
one containing low-gravity, high-luminosity stars — the 
PNNVs RXJ2117-f3412, HE1429-1209, Kl-16, NGC 246, 
HS2324-H3944, Lo 4 and AbcU 43 — , and the other cor- 
responding to high-gravity stars — the naked GW Vir 
stars PG 1159-035, PG 2131+066, PG 1707+427 and PG 
0122+200. Finally, we note that among the low-gravity 
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Fig. 1. The observed location of PG 1159 stars in the logTcs — \ogg plane. Filled squares represent GW Vir pul- 
sators, hollow squares represent non-pulsators PG 1159 stars, diamonds correspond to "hybrid-PG 1159" stars and 
triangles represent stars with no photometric data available. Data are taken from Kawaler et al. (2004), Werner 
et al. (1997), Dreizler & Heber (1998), Nagel & Werner (2004), and Vauclair et al. (2005). PG 1159 evolutionary 
tracks — this work; see also Althaus et al. (2005) — are indicated with solid curves; stellar masses are from right 
to left 0.53, 0.54, . . . , O.64M0. Dashed curve corresponds to an AGB-manque evolutionary track of a model star with 
M, — O.49M0 (Scoccola et al. 2005). In this and the remainig figures displaying results from our additional evolu- 
tionary sequences with different stellar masses, only results for evolutionary stages after the transitory phase has long 
vanished are shown (see text toward the end of 



GW Vir stars there are two objects labeled as pulsating 
"hybrid-PG 1159" — HS2324-h3944 and Abell 43. They 
are variable PG 1159 stars showing strong Balmer lines of 
hydrogen in their spectra. 

A longstanding problem associated with GW Vir stars 
is related to the excitation mechanism for the pulsations 
and, in particular, to the chemical composition of the driv- 
ing zone, an issue that has always attracted the attention 
of the theorists but that until very recently has eluded all 
the attempts of understanding. In their pioneering works, 
Starrfield et al. (1983, 1984, 1985) and StanghcUini et 
al. (1991) realized that g-mode pulsations with periods 
around 500 s could be driven in PG 1159 models by the 
K-mechanism associated with the partial ionization of car- 



bon and oxygen. While these authors were successful in 
finding the correct destabilizing mechanism, the periods 
of the excited modes were too short as compared with the 
observed periodicities in GW Vir stars. In addition, their 
models required a driving region very poor in helium in 
order to be capable to excite pulsations. The latter re- 
quirement led to the suspicion that a composition gradi- 
ent would exist to make compatible the helium-devoid en- 
velope and the helium-rich atmospheric layers. Even very 
recent theoretical work by Bradley & Dziembowski (1996) 
and Cox (2003) points out the necessity of a chemical com- 
position in the driving region different from the observed 
surface composition. This is difficult to explain in view of 
the fact that PG 1159 stars are still experiencing mass loss. 
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a fact that prevents the action of gravitational setthng of 
carbon and oxygen. Clearly at odds with conclusions of 
these works, the calculations by Saio (1996) and Gautschy 
(1997), and more recently by Quirion et al. (2004) and 
Gautschy et al. (2005), demonstrate that ^r-mode pulsa- 
tions in the range of periods of GW Vir stars could be 
easily driven in PG 1159 models at satisfactory effective 
temperatures with an uniform envelope composition com- 
patible with the observed photospheric abundances. 

While non-adiabatic stability issues like those men- 
tioned above are of utmost importance in elucidating the 
cause of GW Vir pulsations and to place constraints on 
the composition of the subphotospheric layers, there is 
other avenue that may help to extract valuable informa- 
tion about the internal structure and evolutionary sta- 
tus of PG 1159 stars, by considering adiabatic pulsa- 
tion periods alone. This approach, known as asteroseis- 
mology, allows us to derive structural parameters of in- 
dividuals pulsators by matching adiabatic periods with 
the observed periodicities. Current examples of astero- 
seismic studies of GW Vir stars are those of Kawaler 
& Bradley (1994) (PG 1159-035), Kawaler et al. (1995) 
(PG 2131-F066), O'Brien et al. (1998) (PG 0122+200), 
Vauclair et al. (2002) (RXJ2117-h341) and Kawaler et al. 
(2004) (PG 1707-F427). These studies have been success- 
ful in providing precise mass determinations and valuable 
constraints on the outer compositional stratification of PG 
1159 stars. As important as they were, these investiga- 
tions were based on stellar models that were not actually 
fully self- consistent evolutionary PG 1159 models. In ad- 
dition, with the notable exception of the exemplary work 
of Kawaler & Bradley (1994), all these studies used the ob- 
served mean period spacing and the asymptotic theory of 
stellar pulsations alone to make asteroscismic inferences; 
that is, no detailed period fitting were carried out in those 
works. In part, the reason for this is that for most of GW 
Vir stars (with exception of PG 1159-035) the observed 
mode density appear to be not high enough for detailed 
asteroscismic analysis. 

It is important to mention that for applications requir- 
ing accurate values of the adiabatic pulsation periods of 
GW Vir stars only full evolutionary models can be used. 
This is particularly true for asteroscismic fits to individual 
stars. In fact, as it has been shown, stellar models repre- 
senting PG 1159 stars should reflect the thermal structure 
of their progenitors because is that structure which pri- 
marily determines the adiabatic period spectrum, at least 
at the high luminosity phases (Kawaler et al. 1985). This 
is of utmost importance concerning the theoretically ex- 
pected rates of the pulsation period change — see, e.g., 
Kawaler et al. (1986). In addition, special care must be 
taken in the numerical simulations of the evolutionary 
stages prior to the AGB phase. Specifically, the model- 
ing process should include extra mixing episodes beyond 
the fully convective core during central helium burning 
(Straniero et al. 2003). In particular, overshoot during the 
central helium burning stage leads to sharp variations of 
the chemical profile at the inner core — see Straniero et al. 



(2003); Althaus et al. (2003). This, in turn, could have a 
strong impact on the mode-trapping properties of the PG 
1159 models. Last but not least, a self-consistent treat- 
ment of the post-AGB evolution is required. In fact, the 
details of the last helium thermal pulse on the early WD 
cooling branch and the subsequent born-again episode de- 
termines, to a great extent, the chemical stratification and 
composition of the outer layers of the hot PG 1159 stars. 
The composition stratification of the stellar envelope, in 
turn, plays an important role in the mode-trapping char- 
acteristics of PG 1159 stars — see Kawaler & Bradley 
(1994). 

All these requirements are completely fulfilled by 
the O.5895-M0 PG 1159 models presented recently by 
Althaus et al. (2005). In fact, these authors have computed 
the evolution of hydrogen-deficient post-AGB WDs tak- 
ing into account a complete and detailed treatment of the 
physical processes that lead to the formation of such stars. 
These models have been recently used for non-adiabatic 
studies of GW Vir stars by Gautschy et al. (2005). The aim 
of this paper is to explore the adiabatic pulsation proper- 
ties of GW Vir stars by employing the full evolutionary 
PG 1159 models of Althaus et al. (2005). In addition, we 
employ these models to make seismic inferences on sev- 
eral GW Vir stars. The paper is organized as follows. In 
Sect. 2, we describe the main characteristics of our numer- 
ical codes and the PG 1159 models employed. The adia- 
batic pulsational properties for our PG 1159 sequences 
are described in detail in §3. Special emphasis is given 
to the mode-trapping properties of our models. Sect. 4 is 
devoted to describing the application of our extensive pe- 
riod computations for asteroscismic inferences of the main 
structural properties of various high-gravity GW Vir stars. 
Finally, in §5 we briefly summarize our main findings and 
draw our conclusions. 

2. Evolutionary code, input physics and stellar 
models for GW Vir stars 

The adiabatic pulsational analysis presented in this work 
is based on artificial PG1159 evolutionary sequences with 
different stellar massess. These sequences have been de- 
rived from a full post-born again evolutionary sequence 
of 0.5895 Mq (Althaus et al. 2005) that has been com- 
puted taking into account the evolutionary history of the 
progenitor star. Specifically, the evolution of an initially 
2.7 Mq stellar model from the zero-age main sequence has 
been followed all the way from the stages of hydrogen and 
helium burning in the core up to the tip of the AGB where 
helimu thermal pulses occur. We mention that during the 
thermally pulsing AGB phase our models experience ap- 
preciable third dredge-up which causes the carbon-oxygen 
ratio to increase from ^ 0.25 to « 0.34 by the time the 
remnant departs from the AGB. Apart from implications 
for the chemical stratification of the post-AGB remnant, 
the occurrence of the third dredge-up mixing event affects 
the degeneracy of the carbon-oxygen core, which in turn 
affects the location of the PG1159 tracks (see Werner & 
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Fig. 2. The surface (thick line) and hehum-burning (thin 
line) luminosities (logL^/L© and logLnc/i©, respec- 
tively) in terms of the effective temperature. The up- 
per panel depicts the situation for our consistent 0.5895- 
Mq sequence. The bottom panel corresponds to the 0.64- 
Mq artificially generated sequence, in which the the 
unphysical evolutionary stages are indicated with dotted 
lines. For details, see text. 



Herwig 2006 for comments). After experiencing 10 ther- 
mal pulses, the progenitor departs from the AGB and 
evolves toward high effective temperatures, where it expe- 
riences a final thermal pidse (hiring the early WD cooling 
phase — a very late thermal pulse and the ensuing born- 
again episode, see Blocker (2001) for a review — where 



most of the residual hydrogen envelope is burnt. After the 
occurrence of a double-loop in the Hertzsprung-Russell 
diagram, the now hydrogen-deficient, quiescent helium- 
burning O.5895-M0remnant evolves at almost constant lu- 
minosity to the domain of PG 1159 stars with a surface 
chemical composition rich in helium, carbon and oxygen: 
('*He,i2C,i60)= (0.306, 0.376, 0.228)^. This is in good 
agreement with surface abundance patterns observed in 
pulsating PG 1159 stars (Dreizler & Hebcr 1998; Werner 
2001). Also, the surface nitrogen abundance (about 0.01 
by mass) predicted by our models is in line with that de- 
tected in pulsating PG 1159 stars (see Dreizler & Heber 
1998). 

We mention that the nuclear network considered in 

the stellar modelling accounts explicitly for 16 chemi- 
cal elements and 34 thermonuclear reaction rates to de- 
scribe hydrogen and helium burning. Abundances changes 
for the elements considered was included by means of 
a time-dependent scheme that simultaneously treats nu- 
clear evolution and mixing processes due to convection 
and overshooting. A treatment of this kind is particularly 
necessary during the extremely short-lived phase of the 
born-again episode, for which the assumption of instanta- 
neous mixing is inadequate. In particular, overshooting is 
treated as a diffusion process (see, e.g., Herwig et al. 1997) 
and has been considered during all evolutionary phases. 
Radiative opacities are those of OPAL (including carbon- 
and oxygen-rich compositions; Iglesias & Rogers 1996), 
complemented, at low temperatures, with the molecular 
opacities from Alexander & Ferguson (1994). 

Because of the considerable load of computing time 
and numerical difficulty issues involved in the generation 
of PG 1159 stellar models that take into accoimt the com- 
plete evolutionary history of the progenitor star, we were 
forced to consider one case for the complete evolution, that 
is, the 0.5895-A'/q sequence. However, to reach a deeper 
understanding of the various trends exhibited by pulsa- 
tions in GW Vir stars, an assessment of the dependence 
of the pulsation properties of PG1159 models with the 
several structural parameters should be highly desirable. 
This statement is particularly true regarding the stellar 
mass of the models. In this paper we consider a set of 
additional sequences of models within a small range of 
stellar masses. In absence of actual self-consistent evolu- 
tionary computations of post-born again PG1159 models 
with different stellar masses, we elect to create several 
evolutionary sequences by employing an artificial proce- 
dure starting from our "seed" model sequence of 0.5895 
Mq (see below) . This approach has been adopted in many 
pulsation studies aimed at exploring the effects of different 
model parameters on the pulsation properties (sec, e.g., 
Kawaler & Bradley 1994), starting from a post- AGB evo- 
lutionary model. At variance with previous calculations. 



^ In particular, the surface abundance of oxygen depends on 
the free parameter / of overshooting (which is a measure of 
the extent of the overshoot region) and the number of thermal 
pulses considered during the AGB phase (see Herwig 2000). 
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in this work we derived additional model sequences with 
stellar masses in the ranges 0.53 — 0.58 Mq and 0.60 — 0.64 
Mq (with a step of 0.01 Mq ) from the post-born again 
O.5895-M0 sequence previously described. To this end, 
we artificially changed the stellar mass^ to the appropri- 
ate values shortly after the end of the born-again episode. 
Although this procedure led to a sequence of unphysical 
stellar models for which the hehum-burning luminosity is 
not consistent with the thermo-mechanical structure of 
the models, we found that the transitory stage vanishes 
at high luminosities, before the star reaches the "knee" 
in the log Toff — logg plane'^ This can be understood by 
examining Fig. [21 which shows the surface and helium- 
burning luminosities in terms of the effective temperature. 
The upper panel depicts the situation for our consistent 
0.5895-Mq sequence. The bottom panel corresponds to 
the O.64-M0 artificially generated sequence. Note that, 
as result of the change in the thermal structure of the 
outer layers induced by the abrupt change in the stellar 
mass (at low effective temperatures), the helium-burning 
luminosity becomes seriously affected. Note that once the 
transitory stage has disappeared, the helium-burning lu- 
minosity follows a similar trend to that expected from 
the full evolutionary computation (upper panel). In the 
bottom panel, the unphysical evolutionary stages are in- 
dicated with dotted lines, while the solid lines correspond 
to the evolutionary stages which we consider valid for the 
purposes of this paper. Thus, the tracks depicted in Fig.^ 
correspond to those stages of the evolutionary sequences 
for which the unphysical transitory has long disappeared. 
We stress that pulsation results presented along this work 
will be performed on stellar models belonging to the evo- 
lutionary sequences shown in Fig. (see also Fig.O. We 
warn again that the pulsational results presented here are 
based on stellar models artificially constructed, and that 
full evolutionary models with different stellar masses and 
complete evolutionary history are required to place such 
results on a solid background. 

We stress that, since we do not account for the com- 
plete evolution of the progenitor stars (except for the 
O.5895-M0 sequence), the internal composition of the 
models we use in part of our pulsational analysis may 
not be completely realistic, a fact that should be kept in 
mind when we analyze, in particular, the dependence of 
the asymptotic period spacing and mode trapping prop- 

^ We change the stellar mass in the following way. In our evo- 
lutionary code, the independent variable is ^ = ln(l — rrir/M,), 
where M, is the stellar mass and rrir is the mass contained in 
a sphere of radius r. When the stellar mass is changed (say to 
Ml), the ^ values at each grid point are the same as before; 
so, the new mass at a given ^ is = Ml{l — exp(^)). The 
chemical composition at each ^ is the same as before, but the 
mass contained at ^ is different. For instance, if Ml < M*, 
then mj, < rrir at a given ^. Note that this procedure is differ- 
ent from, for instance, simply extracting mass from the outer 
layers. 

^ This procedure has been employed, for instance, by O'Brien 
(2000) on a O.573-M0post-AGB evolutionary model. 
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Fig. 3. Evolutionary tracks for sequences with stellar 
masses, from right to left, of 0.53, 0.54, . . . , O.64M0. The 
sequence with O.5895M0 is indicated with thicker line. 
For this sequence, the black dots mark the location of five 
specific models to be analyzed in ti3.2l 

erties with the stellar mass. We also stress that the shape 
of the carbon/oxygen profile — particularly in the core 
— is quite sensitive to the adopted reaction rates during 
the central helium burning phase (see Kunz et al. 2002; 
Straniero et al. 2003) and the efficiency of the convection- 
induced mixing (overshooting and/or semiconvection; see 
Straniero et al. 2003). In particular, we have recently 
found (Corsico & Althaus 2005) that with a lower rate 
of the ^^C(a,7)^^0 reaction than we adopted here, the 
central abundances of oxygen and carbon become quite 
similar, partially smoothing the chemical steps in the core. 
However, as shown in that paper, the mode trapping struc- 
ture resulting from the core is not seriously altered. 

3. Adiabatic pulsational analysis 

In this section we present an extensive adiabatic analysis 
of our grid of PG 1159 evolutionary models. Since the most 
detailed and thorough reference in the literature on adia- 
batic pulsations of PG 1159 stars is the work of Kawaler & 
Bradley (1994) [hereinafter KB94], we shall invoke repeat- 
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edly that paper throughout this work to compare their 
results with our own findings. We have assessed the pul- 
sation properties of about 3000 stellar models, the stellar 
mass of which range from 0.53 Mq to 0.64 Mq, embrac- 
ing the well-established mean mass for WDs. Fig.|21shows 
the evolutionary tracks of our model sequences. All of our 
PG 1159 sequences are characterized by a helium layer 
thickness of O.OO52M0 — as predicted by our evolution- 
ary calculations for the O.5895-Af0 models (Althaus et al. 
2005). For each analyzed model we have computed (7- mode 
periods in the range 50 — 2500 s, although in some cases 
we have extended the range of periods and also included 
p-mode computations. We limited our calculations to the 
degrees £=1,2 because the periods observed in pulsating 
PG 1159 stars have been positively identified with ^=1,2 
(see, e.g., Winget et al 1991; Vauclair et al. 2002). Since 
one of the purposes of this work is to get values of adi- 
abatic periods as precise as possible, we have employed 
about 2700 — 3000 mesh-points to describe our background 
stellar models. 

3.1. Pulsation code and adiabatic quantities 

We begin this section with a brief description of the pulsa- 
tion code employed and the relevant pulsation quantities. 
For computing pulsation modes of our PG 1159 models we 
use an updated version of the pulsational code employed 
in Corsico et al. (2001). Briefly, the code, which is coupled 
to the LPCODE evolutionary code, is based on the gen- 
eral Newton- Raphson technique and solves the full fourth- 
order set of equations governing linear, adiabatic, nonra- 
dial stellar pulsations following the dimensionless formula- 
tion of Dziembowski (1971). The pulsation code provides 
the dimensionless eigenfrequency Uk — k being the ra- 
dial order of the mode — and eigenfunctions yi, • • ■ , 2/4. 
From these basic quantities the code computes the pulsa- 
tion periods (11^), the oscillation kinetic energy (Kk), ro- 
tation splitting coefficients (Cfc), weight functions {Wk), 
and variational periods (Hp for each computed eigen- 
mode. Usually, the relative difference between 11^ and Hk 
is lower than « 10^"*. The pulsation equations, boundary 
conditions and relevant adiabatic quantities are given in 
Appendix A. 

The prescription we follow to assess the run of the 
Brunt- Vaisala frequency (N) is the so-called "Ledoux 
Modified" treatment — see Tassoul et al. (1990) — appro- 
priately generalized to include the effects of having three 
nuclear species (oxygen, carbon and helium) varying in 
abundance (see KB94). In this numerical treatment the 
contribution to N from any change in composition is al- 
most completely contained in the Ledoux term B; this 
fact renders the method particularly useful to infer the 
relative weight that each chemical transition region have 
on the mode-trapping properties of the model (see M3.4() . 
Specifically, the Brunt- Vaisala frequency is computed as: 

^2 XT 
P Xp 



where 

1 dl 

XT^ d 



being 





d\np 




d\np 




d\np 


XT = 


d\nT 


Xp = 

p 




XX. = 

T 





(3) 



Fig.0]shows a representative spatial run of the Brunt- 
Vaisala frequency for a PG 1159 model. The model is 
characterized by a stellar mass of 0.5895 Mq, an effec- 
tive temperature of « 139 000 K and a luminosity of 
log(L*/LQ) = 2.31. In addition, the plot shows the in- 
ternal chemical stratification of the model for the main 
nuclear species (upper region of the plot) and for illus- 
trative purposes the profile of the Ledoux term B (inset). 
The figure emphasizes the role of the chemical interfaces 
on the shape of the Brunt- Vaisala frequency. In fact, each 
chemical transition region produces clear and distinctive 
features in N , which are eventually responsible for the 
mode-trapping properties of the model. At the core re- 
gion there are several peaks at M^jM^ w 0.4 — 0.6 re- 
sulting from steep variations in the inner oxygen/carbon 
profile. The stepped shape of the carbon and oxygen abun- 
dance distribution within the core is typical for situations 
in which extra mixing episodes beyond the fully convective 
core during central helium burning are allowed. In partic- 
ular, the sharp variation around My w 0.56M*, which is 
left by overshoot at the boundary of the convective core 
during the central helium burning phase, could be a poten- 
tial source of mode-trapping in the core region — "core- 
trapped" modes; see Althaus et al. (2003) in the context 
of massive DA WD models. The extended bump in N at 
My. K, 0.96M, is other possible source of mode-trapping, 
in this case associated with modes trapped in the outer 
layers. This feature is caused by the chemical transition of 
helium, carbon and oxygen resulting from nuclear process- 
ing in prior AGB and post-AGB stages. In a subsequent 
Section we shall explore in detail the role played by the 
chemical transition regions on the mode-trapping proper- 
ties of our PG 1159 evolutionary models. 

3.2. The basic nature of the nonradial pulsation 
spectrum 

This section is devoted to exploring at some length the 
basic properties of nonradial pulsation modes of our PG 
1159 pre- WD evolutionary models. To this end, we have 
extended the scope of our pulsation calculations by includ- 
ing low order g-modes and also nonradial p-modes in the 
discussion, although we are aware that short periods like 
those associated with such modes have not been ever de- 
tected in variable PG 1159 stars. Here, we have chosen to 
analyze the sequence of O.5895M0-models from the early 
phases of evolution at constant luminosity — shortly be- 
fore the maximum of effective temperature — until the 
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M/M. 

Fig. 4. The run of the squared Brunt- Vaisala frequency in 
terms of the mass coordinate, corresponding to a 0.5895- 
MqFG 1159 model at Toff = 139 000 K and log(L*/L0) = 
2.31. The mass of the hehum content is of 0.008821M* 
(0.0052Mq). Dark regions denote the contributions of the 
Ledoux term B (shown in the inset) to the Brunt- Vaisala 
frequency. The chemical profile is displayed in the upper 
zone of the plot. 



beginning of the WD cooling track — when models evolve 
at roughly constant radius. We stress that this is the first 
time that the adiabatic pulsation properties of fully evo- 
lutionary post born-again PG 1159 models are assessed 
— an exception is the non- adiabatic analysis of Gautschy 
et al. (2005). Thus, we believe that it is instructive to 
look into the basic features of the full nonradial pulsation 
spectrum of these models. 

We begin by describing the evolution of the period 
spectrum. For the O.5895M0 sequence we have available 
physically meaningful stellar models from stages previous 
to the beginning of the WD cooling track. This fact al- 
lows us to perform an exploratory exercise of the period 
spectrum from the moment when the objects still retain 
the thermal structure resulting from prior evolution of 
the progenitor star. In Fig. we show the £ = 1 pul- 
sation periods 11^ in terms of the effective temperature 
before (panel a) and after (panel b) the model reach the 
turning point in the H-R diagram (see Fig. We find 
the pulsation periods to decrease in stages in which the 
model is still approaching to its maximum effective tem- 
perature (Tcff « 176 600 K). This is because the star is 
undergoing a rapid contraction, in particular at the outer 
layers. As is well known, contraction effects lead to a de- 
crease in the pulsation periods (see Winget et al. 1983). 
When the model has already settled into their cooling 
phase, the periods increase in response to the decrease 
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Fig. 5. The evolution of the £ = 1 pulsation periods 11^ in 
terms of the effective temperature, corresponding to the 
sequence of 0.5895Mo-models. Panel (a) depicts the situ- 
ation before the models reach the maximum value of ef- 
fective temperature (Max. Teff), and panel (b) shows the 
situation when the models have already passed through 
this stage. Arrows indicate the effective temperature cor- 
responding to two elected models that are analyzed in 
Figs. El and Note that the orientation of the T^s axis 
in panel (a) is opposite to that of panel (b). See text for 
details. 

of the Brunt- Vaisala frequency. Thus, the behavior exhib- 
ited is typical of WD pulsators, with increasing periods as 
the effective temperature decreases. Note in panel (a) of 
Fig. Elthat when Tcs < 176 600 K, the low order periods 
(life <^ 200 s) exhibit a behavior reminiscent to the well 
known "avoided crossing" . When a pair of modes experi- 
ences avoided crossing, the modes exchange their intrin- 
sic properties (see Aizenman et al. 1977). As we shall see 
below, these are modes with a mixed character, that is, 
modes that behave as g-modes in certain zones of the star 
and as p-modes in other regions. 

In order to obtain valuable qualitative information of 
nonradial pulsations we apply a local treatment to the 
pulsation equations (Unno et al. 1989). Employing the 
Cowling approximation, in which the perturbation of the 
gravitational potential is neglected (<&' = 0), and assum- 
ing that the coefficients of the pulsation equations are al- 
most constant — that is, in the limit of high radial order 
k — we obtain simple expressions for the eigenfunctions: 
yi,?/2 oc e'}cp{ikrr), and an useful dispersion relation fol- 
lows: k'^ = a~'^c~'^ (cr^ — iV^) (cr^ — Lj), which relates the 
local radial wave number kr to the pulsation frequency a. 
Li is the Lamb frequency, defined as Lij = £{£ + l)c^/r^. 
Note that if > N'^,L'j or a'^ < N'^^L'j, the wave num- 
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Fig. 6. Propagation diagram corresponding to a 
PG 1159 evolutionary model characterized by 
[M,/Mo,log(L*/LQ), log Toff] = [0.5895,3.68,5.168] 
for £ = I modes. The effective temperature of this model 
is marked with an arrow in panel (a) of Fig. |31 White 
(gray) areas correspond to propagation (evanescent) 
zones according the local analysis described in the text. 
The square of the Brunt- Vaisala frequency (N) and the 
Lamb frequency (Lg) are depicted with dashed and solid 
lines, respectively. Horizontal thin lines show the square 
of the eigenfrequencies, and small symbols show the loci 
of the nodes of the radial eigenfunction, yi (filled circles: 
g-modes; filled squares: p-modes). Open circles correspond 
to nodes of modes with mixed character. These modes 
propagate in both P- and G-regions. Note the excellent 
agreement between the full numerical solution and the 
predictions of the simple local treatment. 

ber kr is real, and if N'^ > a'^ > or N'^ < a'^ < L^, kr 
is purely imaginary. 

In Fig.lHlwe show and iV^ as functions of — log(l — 
Mr/m^:). This type of diagrams are called propagation di- 
agrams (see, e.g., Unno et al. 1989; Cox 1980). The figure 
corresponds to a PG 1159 model with A/* = O.5895M0, 
log(L*/L0) = 3.68 and T^b = 147 200 K. This model cor- 
responds to a stage in which the star is still approaching to 
its maximum Toff (see the arrow in panel (a) of Fig. [S]). The 
dashed curve shows the square of the Brunt- Vaisala fre- 
quency, and the solid curve depicts the run of the square 
of the Lamb frequency for = 1. Since this model still 
retains some similarities to those of the red giants from 
which these models are descendants, the Brunt- Vaisala 
frequency is characterized by large central values as com- 
pared to those characterizing the outer regions. In the plot, 
gray areas correspond to regions in which kr is purely 
imaginary. As a result, modes should be unable to prop- 




2 4 6 8 10 



-log(1-M/M.) 

Fig. 7. Same of Fig. but for a PG 1159 model 
characterized by [Af,/M0, log(L,/i0), log Toff] = 
[0.5895,2.31,5.143]. The effective temperature of the 
model is that marked with an arrow in panel (b) of Fig. 
|S1 Note that in this case the G and P branches of modes 
are clearly separated in frequency. 



agate in these "evanescent" zones. In contrast, in regions 
in which kr is real (white areas in the plot) modes should 
be oscillatory in space, and these regions correspond to 
"propagation" zones. As we can see in the figure, there 
exist two propagation regions, one corresponding to the 
case (7^ > N^, Lj, associated with p- modes (P-region), and 
other in which the eigenfrequencies satisfies < N'^, Lf, 
associated with (/-modes (G- region). The predictions of 
this simple local analysis are nicely confirmed by our full 
numerical solution. Indeed, the figure shows the computed 
eigenfrecuencies of p- and g-modes (thin horizontal lines) 
and the loci of the nodes of the radial eigenfunction, yi 
(small symbols). Note that no node lies in the evanescent 
regions, meaning that the radial eigenfunction is not os- 
cillatory here. This figure clearly shows that the p-modes 
and g-modes are trapped in the P-region and the G-region, 
respectively. 

Note that, because the very peculiar shape of the 
Brunt- Vaisala frequency in the inner regions of the star, 
there is a considerable range of intermediate frequencies 
in which the inner parts of the star lie predominantly in 
the G-region, and the outer parts lie in the P-region. Thus, 
these intermediate- frequency modes possess a mixed char- 
acter: they behave like g-modes in the inner parts of the 
star and like p-modes in the outer parts. In Fig. |S1 these 
modes are indicated with open circles. Note that nodes 
in yi for such modes may occur in both the P- and G- 
regions. In order to get the radial order of these modes 
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Fig. 8. The normalized weight function of the mode g^o 
vs the stellar radius corresponding to O.5895Af0-models at 
different effective temperatures. The location of the stellar 
models in the H-R diagram is shown in Fig. |3| Note that 
the regions in which Wk adopts appreciable amplitude — 
the regions of period formation — move from the core 
toward the surface of the model. 

we employ the scheme of Scuflaire (1974)'^. In the figure 
modes with mixed character and some low order "pure" 
p- and g-modes are appropriately labeled. 

Fig. [71 shows a propagation diagram corresponding to 
a model characterized by Af, = O.5895M0, log(L,/Lo) — 
2.31 and Toff = 139 000 K. Unlike the model described in 
Fig- El this PG 1159 model is already located at the early 
WD cooling track — its effective temperature is marked 
with an arrow in panel (b) of Fig. El We expect that the 
fundamental properties of its pulsational spectrum should 
be reminiscent to that of the pulsating WD stars. In fact, 



Specifically, a mode is classified as g if the number of nodes 
laying in the G-region exceeds the number of nodes laying in 
the P- region (nc > np), and is classified as p if no < np. 
We then assign the radial order according to fc = nc — np for 
(7-modes, or k = np — no for p-modes. 



the Brunt- Vaisala frequency shows larger values in the 
outer regions than in the core, in such a way that the 
propagation regions P and G are markedly separated. As 
a result, modes with mixed character no longer exist, and 
the families of g- and p-modes are clearly differentiated. In 
this case, the radial order of the modes is simply the num- 
ber of the nodes of the radial eigenfunction. Note that for 
this model the g-modes are free to propagate throughout 
the star, at variance with the situation of the model de- 
scribed in Fig. El in which g-modes are strongly confined 
to the core region. 

The regions of the star that most contribute to the pe- 
riod formation are inferred from Fig. El which shows the 
normalized weight function (Eq. (|A.14() of the Appendix) 
of a 5-mode with k = 40, £ = 1 in terms of the stellar 
radius corresponding to 0.5895MQ-models at different ef- 
fective temperatures. The location of these models in the 
H-R diagram is depicted as black dots in Fig.El Note that, 
for the model at TcS « 147 000 K (panel a), Wk adopts 
appreciable amplitudes only in central regions. Since the 
behavior of the weight function for this mode is represen- 
tative of all g-modes with intermediate and high radial 
order fc, we conclude that the g-mode pulsation periods 
for PG 1159 models that have not still reached the maxi- 
mum in TcS are mostly determined by the core regions of 
the star, irrespective of the structure of the outer layers. 
Thus, for models at high-luminosity phases the adiabatic 
pulsations properties reflect the conditions in the degen- 
erate core. The situation for a model placed exactly at the 
turning point in Toft is depicted in panel (b) of Fig.El cor- 
responding to Toff = 176 500 K and log(L*/L0) = 3.26. 
Note that, at this stage, the regions of period formation 
have extended towards more external zones. By the time 
the model has settled into its WD cooling track, the weight 
function exhibits appreciable values throughout the star, 
having some amplitude also in the outer regions, as dis- 
played in panel (c). Below ~ 100 000 K, the weight func- 
tion adopts its maximum value in the outer layers, as can 
be seen in panels (d) and (e). At these stages, the adiabatic 
periods are mainly determined by the outer envelope, al- 
though the inner regions of the star also contribute to 
establishing the period of oscillation. Our findings are in 
full agreement with the results of the early computations 
of Kawaler et al. (1985) — for instance, compare our Fig. 
El with their figure 5. 

3.3. Asymptotic period spacing 

In the asymptotic limit of very high radial order k (long 
periods) , the periods of y-modes of a chemically homoge- 
neous stellar model for a given degree £ and consecutive 
k are separated by a constant period spacing Afl^, given 
by (Tassoul et al. 1990): 



Am = 



Ho 



27r2 



N 
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Fig. 9. The evolution of Hq in terms of stellar lumi- 
nosity for evolutionary models of different stellar masses 
(0.53, ■ • ■ , O.64M0 from top to bottom). As in Figs. [T] and 
only results corresponding to models whose structure 
has completely forgotten the unphysical transitory stage 
(after we switch the stellar mass value; see toward the end 
of Sect. 121) are shown. 

Note that XIq (and An|) is a function of the struc- 
tural properties of the star through the Brunt- Vaisala fre- 
quency. The integral is taken over the g-mode propaga- 
tion region. While for chemically homogeneous models the 
asymptotic formula constitutes a very precise descrip- 
tion of their pulsational properties, the computed g-mode 
period spacings (Allfe = Hk+i — Hfc) in chemically strat- 
ified PG 1159 models show appreciable departures from 
uniformity caused by mode trapping (see next Section). 
However, the average of the computed period spacings 
Allfc is very close to the asymptotic period spacing An| 
when averaged over several modes. 

The evolution of Hq vs the stellar luminosity for mod- 
els with different stellar masses is displayed in Fig. |51 In 
line with previous works (see, e.g., KB94), our results in- 
dicate that Ho increases with decreasing stellar mass. In 
particular, the value of Hq exhibits an increment of about 
6 — 8 s when the stellar mass decreases about O.IA/q. This 
is due to the fact that the asymptotic period spacing is in- 
versely proportional to the integral of TV (Eq. 01 . As we 
go to lower stellar masses, the Brunt- Vaisala frequency 
decreases and consequently IIq goes up. 

Fig. shows that IIq is also a function of the stellar 
luminosity, although the dependence is much weaker than 
the dependence on the stellar mass. As documented by 
the figure, for \og{L^/LQ) <^ 3.25 — 2.75 — depending 
on the stellar mass — the period spacing increases with 
decreasing the luminosity. This fact can be explained on 



the basis that, due to the increasing degeneracy in the core 
as the star cools, the Brunt- Vaisala frequency gradually 
decreases, causing a slow increment in the magnitude of 
Hq . Finally, since in this work we do not consider different 
helium layer mass values, we have not been able to assess 
the possible dependence of the asymptotic period spacing 
on Mjie- However, this dependence is expected to be very 
weak (see KB94). 

3.4. Mode-trapping properties 

As we have mentioned before, the period spectrum of 
chemically homogeneous stellar models is characterized by 
a constant period separation, given very accurately by Eq. 
0). However, current evolutionary calculations predicts 
that PG 1159 stars — as well as (hydrogen-rich) DA and 
(helium-rich) DB WD stars — must have composition gra- 
dients in their interiors, something that also the observa- 
tion indicates. The presence of one or more narrow regions 
in which the abundances of nuclear species are rapidly 
varying strongly modifies the character of the resonant 
cavity in which modes should propagate as standing waves 
— the propagation region. Specifically, chemical interfaces 
act like reflecting boundaries that partially trap certain 
modes, forcing them to oscillate with larger amplitudes in 
specific regions — bounded either by two interfaces or by 
one interface and the stellar centre or surface — and with 
smaller amplitudes outside of that regions. The require- 
ment for a mode to be trapped is that the wavelength of 
its radial eigenfunction matches with the spatial separa- 
tion between two interfaces or between one interface and 
the stellar centre or surface. This mechanical resonance, 
known as mode trapping, has been the subject of intense 
study in the context of stratified DA and DB WD pul- 
sations — see, e.g.. Brassard et al. (1992); Bradley et al. 
(1993); Corsico et al. (2002). In the field of PG 1159 stars, 
mode trapping has been extensively explored by KB94; we 
refer the reader to that work for details. 

There are (at least) two ways to identify trapped 
modes in a given stellar model. First, we can consider the 
oscillation kinetic energy, K^. We stress that, because the 
amplitude of eigenfunctions is arbitrarily normalized at 
the model surface in our calculations (see the Appendix), 
the values of the kinetic energy are useful only in a rela- 
tive sense. Kk is proportional to the integral of the squared 
amplitude of eigenfunctions, weighted by the density (see 
Eg. IA.131 in the Appendix). Thus, modes propagating in 
the deep interior of the star, where densities are very high, 
will exhibit larger values than modes which are oscillating 
in the low-density, external regions. When only a single 
chemical interface is present, modes can be classified as 
modes trapped in the outer layers, modes confined in the 
core regions, or simply "normal modes" — which oscil- 
late everywhere in the star — characterized by low, high, 
and intermediate values, respectively (see Brassard et 
al. 1992). This rather simple picture becomes markedly 
complex when the stellar model is characterized by sev- 
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Fig. 10. The forward period spacing Allfc vs period 11^ for 
dipole = 1) modes. The asymptotic period spacing An| 
as given by Eq. I^J is depicted with a thin horizontal line. 
Panel (a) corresponds to the same PG 1159 model consid- 
ered in Fig. ^ Remainder panels correspond to the same 
stellar model but for the situation in which the Ledoux 
term has been artificially suppressed in specific regions of 
the model: (b) B = at M^/M* « 0.96, {c) B = at 
Mr/M^ « 0.4-0.6, and {d) B = everywhere in the star. 



eral chemical composition gradients — see Corsico et al. 
(2002) for the case of ZZ Ceti models. 

A second — and more important from an observational 
point of view — consequence of mode trapping is that the 
period spacing Allfc, when plotted in terms of the pul- 
sation period 11^;, exhibits strong departures from unifor- 
mity. It is the period difference between an observed mode 
and adjacent modes (Afc — ±1) that is the observational 
diagnostic of mode trapping — at variance with Kk, whose 
value is very difficult to estimate only from observation. 
For stellar models characterized by a single chemical in- 
terface, local minima in Allfc usually correspond directly 
to modes trapped in the outer layers, whereas local max- 



ima in Allfc are associated to modes trapped in the core 
region. 

We now turn out to the specific case of our PG 1159 
pre- WD models. As Fig. 01 shows, they exhibit several 
chemical interfaces, some associated with the various steps 
in the ^^O and ^^C profile at the core region, and a single, 
more external transition region in which oxygen, carbon 
and helium are continuously varying (see Fig.^. As stated 
in Sect. 13.11 these composition gradients produce pro- 
nounced peaks in the Brunt- Vaisala frequency — through 
the Ledoux term — that strongly disturb the structure of 
the period spectrum. 

The influence of the chemical composition gradients in- 
side our PG 1159 models on their period spectra is clearly 
shown in panel (a) of Fig. EH in which the £ — 1 pe- 
riod spacing is plotted in terms of the periods for a refer- 
ence model characterized by a stellar mass of O.5895M0, 
an effective temperature of 139 000 K and a luminosity 
of 204. 2Lo. The plot shows very rapid variations of Allfc 
everywhere in the period spectrum, with "trapping ampli- 
tudes" up to about 6 s and an asymptotic period spacing 
of « 20.94 s. The rather complex period-spacing diagram 
shown by Fig. ^| is typical of models characterized by 
several chemical interfaces. In order to isolate the effect of 
each chemical composition gradient on mode trapping, we 
follow the procedure of Charpinet et al. (2000) for models 
of sdB stars — see also Brassard et al. (1992) for the case 
of ZZ Ceti stars. Specifically, we minimize — although no 
completely eliminate — the effects of a given chemical in- 
terface simply by forcing the Ledoux term B to be zero 
in the specific region of the star in which such interface is 
located. In this way, the resulting mode trapping will be 
only due to the remainder chemical interfaces. In the inter- 
est of clarity, we label as "Bochc" the contribution to B 
due to the 0/C/He chemical interface at Mr/M^, « 0.96, 
and "Soc" the contribution to B associated with the 
0/C chemical interface at M^/M* w 0.4 — 0.6 — see 
inset of Fig. ^ Specifically, we have recomputed the en- 
tire g-mode period spectrum under the following assump- 
tions: (1) -BocHo = and Bqc 0- (2) -Bqchc ^ and 
Boc = 0, and (3) Bqcho = and Bqc = (B = in 
all regions). The results are shown in panels (b), (c) and 
(d) of Fig. ^1 respectively. By comparing the different 
cases illustrated, an important conclusion emerges from 
this figure: the chemical transition region at Mr « 0.96M* 
is responsible for the non-uniformities in Allfc only for 
Ilfc <; 500 s (panel c), whereas the chemical composition 
gradients at the core region {Mr « 0.4 — 0.6M*) cause the 
mode-trapping structure in the rest of the period spec- 
trum (panel b) . When B = everywhere inside the model 
(panel d), the period-spacing diagram is characterized by 
the absence of strong features of mode trapping, although 
some structure remains, in particular for low order modes. 
Note the nice agreement between the numerical computa- 
tions of Allfc and the predictions of the asymptotic theory 
given by Eq. 10} for the hmit fc ^ 1. 

From the above discussion, it is clear that the mode- 
trapping features characterizing our PG 1159 models 



12 



Corsico & Althaus: The GW Vir variable stars 




I I I I I J : I I I I I I I 

200 400 600 800 1000 1200 1400 



Fig. 11. Period-spacing diagram and kinetic energy distri- 
bution (upper and lower panel, respectively) for the same 
PG 1159 model considered in panel (a) of Fig.^| Vertical 
dotted lines link local minima in Allfc with the corre- 
sponding values of K^. Maximum values of kinetic energy 
are indicated with arrows. Also labeled are the modes with 
A; = 56 and k = 58, to be analyzed in Fig.lT^ 



are inflicted by the stepped shape of the carbon/ oxygen 
chemical profile at the core — left by prior convective 
overshooting — at least for the range of periods observed in 
GW Vir stars. On their hand, the more external chemical 
transition has a minor influence, except for the regime of 
short periods. We mention that this situation is more ev- 
ident for the more massive PG 1159 models than for the 
less massive ones. 

The findings outlined above are clearly at odd with 
previous results reported by KB94. Indeed, these authors 
have found that the mode-trapping properties of their PG 
1159 models are fixed mainly by the outer 0/C/He tran- 
sition region, to a such a degree that they have been able 
to employ mode-trapping signatures as a sensitive locator 
of this transition region. The origin of this discrepancy 
can be found in the details of the input physics employed 
in build up the background stellar models. Of particular 
interest here is the presence of a much less pronounced 
chemical transition in the outer parts of the C/0 core 
of the KB94 models, as compared with the rather abrupt 
chemical gradients at ~ 0.4 — 0.6Af* characterizing our 
PG 1159 models. If we artificially minimize the effect of 
this transition region — by setting -Bqc = — we imme- 
diately recover the results of KB94. Indeed, the period- 
spacing diagram of panel (c) in our Fig. ^| looks very 
similar to that shown in Fig. 3 of KB94, corresponding to 
Mne/M* = 0.00622. We note, however, that the "trapping 
cycle" — the period interval between two period-spacing 



minima — of our modified model (c) is of « 70 s, whereas 
the value of the KB94 model is « 95 s. This difference is 
due mainly to the fact that our model has a more massive 
hehum envelope (Mho/M* = 0.00881) than that of KB94. 

By means of a simple numerical experiment we have 
identified the main source of mode trapping as due to the 
step-like chemical transition region located at the core. We 
now return to panel (a) of Fig. ^1 and note that there is 
a kind of "beating" which modulates the trapping ampli- 
tude. We also note that this striking feature (seen in all of 
our models) persists even in the case in which the effect 
of the 0/C/Hc transition is artificially minimized, as it 
is illustrated in panel (b). We have found that the beat- 
ing exhibited by the period-spacing distribution is due to 
the combined mode-trapping effects caused by the various 
steps in the 0/C chemical profile in the core. In fact, by 
performing period computations in which only the largest 
peak of the Ledoux term at Mr « 0.58M* (see inset of 
Fig.0J is considered, the beating effect virtually vanishes, 
and the trapping amplitude becomes nearly constant. 

Finally, to understand in more detail the mode- 
trapping properties of our models, we resort to the kinetic 
energy of modes. In Fig.^Jwe plot the period spacing and 
the kinetic energy distribution (upper and lower panel, 
respectively) for the same PG 1159 model considered in 
panel (a) of Fig. ^1 We connect with dotted lines each 
minimum in Allfe with their corresponding value of Kk- 
We note that in spite of the strong non-uniformities ex- 
hibited by AII/j, the Kk values of adjacent modes seem 
to be not very different amongst them. However, a closer 
inspection of the figure reveals that in most of the cases 
a minimum in Allfc corresponds to a local maximum in 
Kk- Our results strongly suggest that modes associated 
with minima in the period- spacing diagram usually corre- 
spond to maxima in the kinetic energy distribution. They 
should be modes characterized by relatively high amplitude 
of their eigenf unctions and weight functions in the high- 
density environment characterizing the stellar core. 

This is precisely what we have found after a careful 
examining of the eigenfunctions at the deepest regions of 
the core. In panels (b) and (c) of Fig.^jwe show the loga- 
rithm of the absolute value of the eigenfunctions yi and 2/2 
(see Eg. lA. Sl ot the Appendix for their definition), respec- 
tively, in terms of the dimcnsionless radius x — r/R^ for 
a mode having a local maxima in kinetic energy (fc = 58, 
thin line) and for a mode having a local minima in Kk 
{k = 56, thick line). Note that below r/Rw 0.3, the eigen- 
functions of the A: = 58 mode have larger amplitudes than 
that of the fc = 56 mode, explaining why the fc = 58 mode 
has a relatively large oscillation kinetic energy. Modes like 
the fc = 58 one are partially trapped in the core region, 
below the O/C chemical interface (see panel a). As stated 
before, these modes correspond generally to minima in the 
period-spacing diagram. In panels (d) and (e) of Fig. [T^ 
the normalized weight function {Wk) is displayed for the 
fc = 56 and fc = 58 modes, respectively. As we mentioned 
previously, the relative values of Wk for a given pulsa- 
tion mode provide information about the specific regions 
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Fig. 12. The spatial run of the modulus of the eigenfunc- 
tions 2/1 and 2/2, shown in panels (b) and (c), respectively, 
and the weight function Wk depicted in panels (d) and 
(e), for dipole pulsation modes with fc = 56 (thick line) 
and fc = 58 (thin line) corresponding to the same PG 1159 
model analyzed in Fig. ^] We also include in panel (a) 
the Ledoux term B and the abundance (by mass) of ^^O, 
12 C and ^He. 



inside the star that most contribute to the period forma- 
tion. From Fig. El we note that for the fc = 58 mode the 
maximum of Wk is identified with the 0/C chemical in- 
terface at r/R « 0.3, and the largest amplitude portion 
of the weight function is located in the deepest regions 
of the core {r/R <^ 0.3). For values of r/R slightly larger 
than 0.3, the weight function abruptly diminishes and then 
reaches a secondary maximum in the region in which ^^C 
becomes more abundant than ^^O. Above this region Wk 
exhibit other secondary maximum at the 0/C/He tran- 
sition zone {r/R « 0.58), and above of this the weight 
function adopts rather lower values {Wk <j 0.3). Thus, 
the shape of the weight function for the A: = 58 mode 
strongly suggests that the deepest regions of the core have 
the larger impact on the period formation of this mode. 
This conclusion is, not surprisingly, in complete agreement 
with the fact that this mode is characterized by a relatively 
high value of the oscillation kinetic energy. The situation 
is markedly different for the case of the A; = 56 mode, 
as it is clearly demonstrated by panel (d) of the figure. 
For this mode, most regions of the star appreciably con- 
tribute to the formation of its period. In particular, the 




Fig. 13. The period-spacing distributions for PG 1159 
models with stellar masses, from bottom to top: 
0.53, 0.54, . . . , O.64M0 at T^s « 140 000 K. In the interests 
of clarity, each curve has been arbitrarily shifted upwards, 
starting from the lowest one, corresponding to the model 
with = O.53M0. Models with stellar masses other than 
0.5895-MQCorrespond to evolutionary stages after the un- 
physical transitory phase (induced by the scaling of the 
stellar mass) has long vanished. 



most important contributions occur in regions located at 
0.2 ^ r/R ^ 0.6. 

3.5. The effect of changes in the stellar mass 

In the above discussion we have been able to identify the 
physical origins and nature of mode trapping in our PG 
1159 models. We wish now to explore the effect of chang- 
ing the stellar mass on the period-spacing diagrams 
and mode trapping. As mentioned previously, we have em- 
ployed in this project model sequences with stellar masses 
in the ranges 0.53 - 0.58 Mq and 0.60 - 0.64 Mq with a 
step of 0.01 Mq , in addition to the 0.5895-Mo sequence 
previously described. We recall that those sequences were 
obtained from the 0.5895-Mq sequence by artificially 
changing the stellar mass according the procedure out- 
lined in [21 Thus, we caution that the results derived in 
this section could change if the complete evolution of pro- 
genitor stars with different stellar masses were accounted 
for. From the asymptotic theory outlined in i l3.3l we know 
that An^ decreases when the stellar mass increases, as it is 
evident from Fig. |^ As mentioned previously, the asymp- 
totic period spacing is very close to the average of the 
computed period spacings Allfc. As a result, the average 
of the computed period spacing also must decreases as the 
stellar mass increases. In fact, Allfc decreases from 23.19 s 
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to 18.69 s when we increase the stehar mass from 0.53Mq 
to O.64M0 for a fixed Toff ~ 140 000 K. 

We also expect variations in the amplitude and cy- 
cle of trapping when we consider changes in the stellar 
mass. In Fig. 1131 we depict period-spacing diagrams for 
PG 1159 models with values of Af* covering the complete 
set of stellar masses considered in this work, and for a 
fixed effective temperature of w 140 000 K. Note that the 
trapping amplitude decreases when the stellar mass in- 
creases. In fact, the maximum of the trapping amplitude 
is of about 10 s for O.53M0, whereas it assumes a value of 
« 5 s, at most, for G.64Mq. This is a direct consequence of 
an increased electron degeneracy in the core of our more 
massive models, which in turn produces a weakening in 
the efficiency of mode trapping. Other striking feature 
observed in Fig. ^| is the decrease in the period of the 
trapping cycle and the shift of mode trapping features to 
lower periods as we go to larger values of the stellar mass. 
Specifically, the portions of the period-spacing diagram 
centered at « 500, « 900, « 1400, « 1700 and « 2200 s for 
the model with M* = 0.53Mq, move to regions centered 
at « 100, ^ 700, 1000, « 1400 and « 1900 s, respec- 
tively, for the model with Af* ~ 0.64^/0. This effect can 
be understood on the basis that, while the 0/C chemi- 
cal structure remains at the same M^/M* value inside the 
model when we consider higher stellar masses, its loca- 
tion in terms of the radial coordinate (r) shifts away from 
the stellar centre, from r/R « 0.22 for = O.53M0 to 
r/R « 0.33 for = O.64M0. A similar effect has been 
found by Corsico et al. (2005) in the context of crystalliz- 
ing ZZ Ceti models, by considering a fixed value of but 
with the inner turning point of the g-mode eigenfunctions 
— the crystallization front — moving outwards. 

3.6. The effect of changes in the effective temperature 

In closing this section, we shall explore the evolution of 
the mode-trapping properties as our PG 1159 evolve to 
cooler temperatures. As we mentioned earlier, the asymp- 
totic period spacing An| and the average of the computed 
period spacings Allfe increase when the stellar luminosity 
decreases (see Fig. EJ. When the star evolves along the 
WD cooling track, a decreasing effective temperature is 
usually associated to a decreasing stellar luminosity. As 
a result, we also expect that An| and Allfc increase as 
cooling proceeds. This trend is confirmed by our pulsa- 
tion computations. In fact, Allfc increases from 20.22 s to 
« 24.12 s when the effective temperature is decreased from 
170 000 K to 70 000 K in a model with O.5895M0. The 
explanation to this effect is straightforward. As the star 
cools their core becomes more degenerate and the Brunt- 
Vaisala frequency decreases because xt goes down (see 
Eq. nj. As a result, the pulsation periods become longer 
and thus the period spacing and the average of the period 
spacing must increase. 

Now, we want to see the effect of changes in the effec- 
tive temperature on the trapping cycle and amplitude. To 
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Fig. 14. Period-spacing distributions for O.5895-Af0 PG 
1159 models at different effective temperatures. In the in- 
terests of clarity, each curve has been shifted upwards, 
except for the lowest one. In a few cases we have included 
numbers at the right side of the plot, corresponding to 
the Tcs value of the models (in units of kK). In partic- 
ular, the thick curve corresponds to the highest value of 
effective temperature reached by the star, T^g ~ 176 kK. 



have a global picture of what happens when the effective 
temperature changes, we depict in Fig. ll4l the evolution of 
the period spacing for O.5895-M0 models, shown at dif- 
ferent effective temperatures (indicated at the right side 
of the plot) starting from a stage in which the model is 
still approaching its maximum Toff (lowest curve) until to 
a phase in which the object is already on the WD cooling 
track (upper curve) . Note that the scale at left of the plot 
has sense only for the lowest curve; the remaining curves 
are shifted upward with purposes of clarity. The figure em- 
phasizes the notable changes that undergoes the period- 
spacing structure as the effective temperature varies. Note 
that there exist two clearly different trends. In stages in 
which the effective temperature increases (before the star 
reaches its hot turning point at Tcff ~ 176, 000 K) the 
mode-trapping features slowly move toward short periods. 
This is due to the fact that pulsation periods decrease in 
response of the overall warming of the star, and that the 
outer layers contract. When the star has already passed 
through its hottest phase, the opposite behavior is exhib- 
ited by the mode-trapping features. Indeed, maxima and 
minima values of the period-spacing distributions move to 
higher periods when the effective temperature diminishes, 
in particular below Toff « 100 000 K. Finally, we note that 
the trapping amplitude increases when the effective tem- 
perature decreases. 
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4. Application to real stars: asteroseismology of 
variable GW Vir stars 

In this section we shall attempt to infer the main struc- 
tural parameters (that is, Af, and Tcff) of several GW 
Vir stars by employing the results of our extensive period 
computations described in previous sections. Due to the 
fact that in this project we have not been able to com- 
pute realistic PG 1159 models with stellar masses other 
than O.5895M0 in the low-gravity phases (before the star 
reaches the "knee" at the highest value of Toff; see Fig. 
nj, we have limited ourselves to make seismic inferences 
only for the high-gravity members of the GW Vir class, 
that is, the naked — stars without surrounding nebulae 

— variables PG 1159-035, PG 2131+066, PG 1707-^427, 
and PG 01224-200. We mention that the results described 
below are conditional on the assumption under the se- 
quences of different stellar masses have been created (see 
Sect.EJ. As it is well known, changes in the thickness of 
the surface helium-rich layer and in the surface composi- 
tion affect the trapping cycle and amplitude (see KB94). 
Thus, by adopting these quantities as free adjustable pa- 
rameters it is possible to perform a fine tunning of the 
trapping cycle and the trapping amplitude to best match 
the period spectrum of a given star. In spite of this, in 
this work we prefer to keep fixed the thickness of the he- 
lium layer and the surface composition by adopting for 
these parameters the values predicted by our evolution- 
ary computations. It is worth to mention that all of these 
naked GW Vir stars exhibit very similar surface abun- 
dances {XB.,Xne,Xc,Xo) « (0.00,0.33,0.50,0.17). Our 
PG 1159 models are characterized by surface abundances 
of (0.000, 0.306, 0.376, 0.228), and thus, they are very 
appropriate for our purposes. Before going to seismic ap- 
plications of our bank of periods, we summarize the basic 
properties of the cited GW Vir stars below. 

The four naked GW Vir stars have higher gravity 
{\ogg J> 7) than the PNNVs, and they are thought to 
be slightly more evolved. From a pulsation point of view, 
the only difference between both types of variable stars 
is that PNNV stars in general pulsate with longer peri- 
ods — in the range of about 1000 — 3000 s — than the 
naked GW Vir — which pulsate with periods below 1000 
s. Specifically, there is a well-defined correlation between 
the luminosity of the star and the periods exhibited: the 
more luminous the star is, the longer the pulsation periods 
are. 

- PG 1159-035: this is the prototype of the PG 1159 spec- 
tral class, and also the prototype of the GW Vir pulsators. 
After the discovery of its variability by McGraw et al. 
(1979), PG 1159-035 became the target of an intense ob- 
servational scrutiny. The most fruitful analysis of its light 
curve was carried out by Winget et al. (1991) employ- 
ing the high-quality data from the Whole Earth Telescope 
(WET; Nather et al. 1990). This analysis showed that PG 
1159-035 pulsate in more than 100 independent modes 
with periods between 300 and 800 s. By using 20 unam- 



biguously identified £ — l^m = periods between 430 
and 817 s, KB94 determined a mass for PG 1159-035 of 
O.59M0, an effective temperature of Toff = 136 000 K and 
a surface gravity of logg = 7.4. In addition, a stellar lumi- 
nosity of 195Lq and a distance of 440 pc from the Earth 
were inferred. On the other hand, spectroscopic analysis 
by Dreizler & Heber (1998) yield T^s = 140 000 K and 
logg = 7, in good agreement with the asteroseismic fit of 
KB94, but a higher luminosity of 5O1L0. 

- PG 2131+066: Was discovered as a variable star by 
Bond et al. (1984) with periods of about 414 and 386 
s, along with some other periodicities. On the basis of 
an augmented set of periods from WET data, Kawaler 
et al. (1995) considered a Tes = 80 000 K and obtained 
a precise mass determination of M^, = O.6IM0, a lumi- 
nosity of IOL0 and a distance from the Earth of 470 pc. 
Spectroscopic constraints of Dreizler & Heber (1998) give 

= O.55M0, Teff = 95 000 K, 39.8L0 and logg = 7.5 
for PG 2131-1-066. By using the updated Dreizler & Heber 
(1998)'s determination of the effective temperature. Reed 
et al. (2000) refined the procedure of Kawaler et al. (1995). 
They found = O.6O8M0, a luminosity of 26L0 and a 
seismic distance to PG 2131-1-066 of 668 pc. 

- PG 1707+427: This star was discovered to be a pulsator 
by Bond et al. (1984). Dreizler & Heber (1998) obtained 
Toff = 85 000 K and logg = 7.5, and a stellar mass and lu- 
minosity of O.54M0 and 25Lq, respectively, were inferred 
from their spectroscopic study. Recently, Kawaler et al. 
(2004) have reported the analysis of multisite observations 
of PG 1707-1-427 obtained with WET. Preliminary model 
fits by using 7 independent £ — 1 modes with periods 
between 334 and 910 ssuggest an asteroseismic mass and 
luminosity of O.57M0 and 23L0, respectively. 

- PG 0122+200: It is the coolest GW Vir variable, with 
Teff = 80 000 K and L* = 20Lq (Dreizler & Heber 1998). 
Spectroscopy indicates a stellar mass of O.53M0 and a 
gravity of logg = 7.5. By employing an analysis based on 
multisite observations with WET, O'Brien et al. (1998) 
report a seismic stellar mass of O.69M0, strikingly higher 
than the spectroscopic estimation. The cooling of PG 
0122-1-200 appears to be dominated by neutrino looses; 
this renders PG 0122-1-200 as the prime target for learn- 
ing neutrino physics (O'Brien et al. 1998). 

In order to infer the structural parameters of these 
four GW Vir stars, we shall employ three different meth- 
ods. First, we estimate the stellar masses by using the 
asymptotic period spacing of our models, as computed 
from Eq. Specifically, we make a direct comparison 
between AH^ and the observed mean period spacing AH, 
assuming that the effective temperature of the target star 
is that obtained by means of spectroscopy. In the second 
approach, we repeat this procedure but this time using the 
average of the computed period spacings AH^ and com- 
paring it with the observed mean period spacing AH. The 
third approach is a fitting method in which we compare 
the theoretical (nj) and observed (Hf) periods by means 
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Table 1. The values of the stellar mass (in Mq) from spectroscopic and asteroseismic studies, and the values obtained 
in this work from three different methods. 

Source PG 1159-035 PG 2131+066 PG 1707+427 PG 0122+200 " 

Spectroscopy 0.54 iO.l'^' 0.55 + 0.1*^' 0.54 + 0.1*^' 0.53 ±0.1^''' 

Asteroseismology 0.59 ± O.Ol'*"' 0.61 ± 0.02''=) 0.57 + 0.02''*) 0.69 ± 0.04^''' 

By comparing All with An? (this work) 0.575 + 0.005 0.615 + 0.015 0.595 + 0.015 0.70 + 0.02 

By comparing ATI with Mh (this work) 0.555 + 0.005 0.575 ± 0.015 0.565 + 0.015 0.64 + 0.02 
By means of period fittings (this work) 0.56 + 0.01 0.60 + 0.01 0.55 + 0.01 0.64 + 0.01 

References; (a) Dreizler & Heber (1998); (b) KB94; (c) Reed et al. (2000); (d) Kawaler et al. (2004); (e) O'Brien (2000) 
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Fig. 15. The asymptotic period spacing for PG 1159 mod- 
els with stellar masses of 0.53, • • • , 0.641/© in terms of 
the effective temperature. This figure is based on Fig. IHl 
The observed mean period spacing AH of the high-gravity 
GW Vir stars PG 1159-035, PG 2131+066, PG 1707+427, 
and PG 0122+200 (21.5 ± 0.1 s, 21.6 ± 0.4 s, 23.0 ± 0.3 s 
and 21.1 ± 0.4 s, respectively) is also shown (taken from 
Kawaler et al. 2004). 

of a standard algorithm. In the three approachs we as- 
sume that the observed periods correspond to £ = l,m = 
modes. 



4.1. Inferences from the asymptotic period spacing, 

An? 

We start by examining Fig. 1151 in which the asymptotic 
period spacing (An|) is shown in terms of the effective 
temperature for PG 1159 sequences with stellar masses of 
0.53, • • • , 0.64Mq. We also include in the figure the values 
of the observed mean period spacing, AH, corresponding 
to the GW Vir stars PG 1159-035, PG 2131+066, PG 
1707+427, and PG 0122+200. The values of ATT and the 
associated error bars are taken from Table 6 of Kawaler 
et al. (2004). In Table 1 we summarize our results. Note 
that the spectroscopic estimation of T^tf for PG 0122+200 
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Fig. 16. The average of the computed period spacings for 
PG 1159 models with different stellar masses in terms of 
the effective temperature. Each panel corresponds to a 
specific GW Vir star. The observed mean period spacing 
An is also shown (taken from Kawaler et al. 2004). 



place it at an effective temperature that is lower than our 
PG 1159 models. Thus, for PG 0122+200 we can only 
make a rough extrapolation. Note the notable agreement 
between our estimates and the values inferred by other 
pulsation studies. Note also that the seismic inferences 
suggest higher mass values as compared with the spec- 
troscopic estimations. We do not actually understand the 
origin of this discrepancy, but we note that spectroscopic 
derivations of the stellar mass are usually very uncertain 
due to the large uncertainty in the determination of log g 
— an error of 0.3 dex translates into an error of O.IM© 
(Dreizler & Heber 1998). Indeed, Dreizler & Heber (1998) 
first estimate logg and T^ff values by using fits to model 
atmospheres, and then they select the stellar mass from 
the evolutionary tracks of O'Brien (2000). Thus, both the 
large uncertainties in the estimation of logg (« 0.5 dex) 
and in the evolutionary computations could account for 
an excessively broad range of allowed masses for a specific 
star. 
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4.2. Inferences from the average of the computed 
period spacings, AUk 

Here, we repeat the procedure described in tl4.1l but this 
time by employing the average of the computed period 
spacings (Allfe) instead the asymptotic period spacing. We 
compare ARk with the measured AH for each star under 
consideration. In computing the average of the computed 
period spacings, we consider only the period interval in 
which the periodicities of a given star are observed. For 
instance, for PG 1159-035, we compute the average of the 
computed period spacings for periods in the range [430, 
841] s. The results of our calculations are shown in Fig ll6l 
in which we show Allfe plotted with dashed lines in terms 
of Tes for the four GW Vir stars under consideration. Note 
that, since we have performed the average of the computed 
period spacings on different period ranges, appropriate to 
the period range exhibited by a specific star, the curves 
of Allfc are different in each panel. In Table 1 we include 
the estimations of the stellar mass for the four stars. Note 
that the stellar masses as determined by this approach 
are appreciably lower than the values inferred by using 
the asymptotic period spacing and thus in better agree- 
ment with spectroscopic inferences. This is due to Allfe 
is typically « 0.8 — 1.0 s smaller than AII^, irrespective 
of the stellar mass and effective temperature. Thus, if we 
compare the observed mean period spacing for a given star 
with Allfe we obtain a smaller total mass (0.02 — 0.06Mq 
lower) than that inferred by comparing the observed mean 
period spacing with AII^. We stress that the asymptotic 
period spacing AII^, as computed by means of Eq. ^ is 
formally valid for the limit of high radial order k in chemi- 
cally homogeneous stars. Because our PG1159 models are 
chemically stratified, we conclude that the estimations of 
from AUk are more realistic than those inferred by 
means of An|. We also found that our values are lower 
than those quoted by other seismic studies (see Table 1). 
We note that, with the exception of the work of KB94 
for PG 1159-035, all these studies compare the observed 
mean period spacing with the asymptotic predictions to 
infer the stellar mass. Thus, not surprisingly, their values 
are very similar to our results obtained from the asymp- 
totic period spacing in ti4.1l but rather departed from our 
values deduced from the average of the computed period 
spacings. 



4.3. Inferences from the computed pulsation periods, 

We turn now to the period fitting procedure. The good- 
ness of a given fit is measured by means of a quality func- 
tion, defined as 



PG 1159-035 



N 



x2(M„Teff)=^min[(np-nJf]. 



(5) 




Here, Hf is one of the N observed pulsation periods, and 
I? 



Hj is a specific computed period with radial index k. 
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T3„[x10^ K] 

Fig. 17. The values as defined by Eq. lO correspond- 
ing to PG 1159-035 versus the effective temperature for 
£ = 1 modes. Each panel depicts the results for a specific 
mass value. The gray strip corresponds to the allowed val- 
ues of Toff according to spectroscopic measurements. The 
arrow indicate the minimum value of associated with 
the effective temperature of the best-fit model. 

The method consists of looking for the PG 1159 model 
that shows the lowest value of x^- Obviously, ^ if 
the match between observed and computed periods were 
perfect. We evaluate x^(Af*,Toff) varying M^^/Mq in the 
range [0.53,0.64] with a step A{M^/Mq) = O.Ol, whereas 
for the effective temperature we adopt a much more finer 
grid, with ATos = 10-20 K. 

We initially searched the optimal model within a 
wide interval of effective temperatures, ranging from the 
hottest point reached by the model in the H-R diagram 
(160 000 - 200 000 K, according to the value of Af,) until 
Tcff ~ 70 000 K. However, throughout our calculations we 
realized that generally the run of x^ as a function of T^s 
for a given star exhibits more than one minimum, meaning 
that the period spectrum of the star could not be fitted 
by an unique PG 1159 model. This behaviour is more pro- 
nounced in the case of the GW Vir stars with few observed 
periods available. This effect can be understood on the ba- 
sis that the pulsation periods of a specific model are gen- 
erally increasing with time. If at a given effective temper- 
ature the model shows a close fit to the observed periods, 
then the function ~x^ reaches a local minima. When the 
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Table 2. Comparison between the observed periods of 
PG 1159-035 (taken from KB94) and theoretical {£ = 1) 
periods of the best-fit model {511 = 1.79 s, asn = 2.21 s). 



PG 2131+066 



1000 







K, 


ATT 


ATT /TT 
Olli/ ill 


100 




Tel 




Tel 


[/o\ 




A QA fl/l 
4oU.U4 


/I QQ 01 


io 


Q 1 7 


n 7/1 
— U. / 4 




4oU. / 1 


A /I o 


iy 


7S 
— Z. / O 


— U.Oz 


10 


4Dy.O ( 


A 7Q c: 1 


on 


— o.y4 


n QA 
— U.o4 




AQA QK 


AGfi KQ 

4yD.0o 


01 

zi 


1 7Q 


— U.oo 




Oi ( .lo 


p;i 7 1 n 


00 
zz 


U.Uo 


n no 
U.Uz 


100 


OOO.iD 


p:q7 c;n 
oo / .OU 


OQ 
ZO 


U.DD 


A 1 O 

U.iz 


ir CO A A 


ooy.lo 


0/1 

z4 


n 71 


— U.lo 




KQi on 


con c^r^ 


oc; 
zO 


n 7/1 
U. / 4 


n 1 Q 
U.lo 


10 


DUO.U4 


oUl.zo 


0(^ 
ZD 


1 7fi 


n oo 
U.zy 


coo fid 


coo cn 
Dzz.DU 


07 
z ( 


u.uu 


u.uu 




fi.A Q /1 1 
D4o.4i 


c/l c 0/1 


OQ 
ZO 


1 QQ 


n OQ 
— U.zc5 


100 


666.22 


665.09 


29 


1.13 


0.17 


687.71 


686.29 


30 


1.42 


0.21 




707.92 


709.95 


31 


-2.03 


-0.29 




729.50 


729.65 


32 


-0.15 


-0.02 


10 


753.12 


750.01 


33 


3.11 


0.41 




773.77 


773.33 


34 


0.44 


0.06 


-^100 


790.94 


795.17 


35 


-4.23 


-0.53 


817.12 


814.96 


36 


2.16 


0.26 




840.02 


836.27 


37 


3.75 


0.45 





model has cooled enough (at some time later), it is possible 
that the accumulated period drift nearly matches the pe- 
riod spacing between adjacent modes (|Afc| — 1). In these 
circumstances, the model is able to fit the star again, as a 
result of which exhibits other local minimum. To break 
this degeneracy, we decided to consider a more restricted 
range of effective temperatures, but that still comfortably 
embraces the spectroscopic estimation of Tcff and its un- 
certainties. Specifically, in searching the "best-fit" model 
for a specific star we considered effective temperatures lay- 
ing inside an interval r'P'='=-20, 000 < T^s < T"P'='=-f 20 000 
K , where T^g""^ is the spectroscopic determination of the 
effective temperature. As we shall describe below, for the 
star PG 1159-035 we were able to find a best-fit model 
for which adopts the lowest value within this interval 
of TcffS. The situation was considerably more difficult for 
PG 2131+066, PG 1707+427, and PG 0122+200, due to 
the scarcity of pulsation periods characterizing their pul- 
sation spectra. In fact, for these stars the ambiguity as- 
sociated with the existence of multiple solutions remained 
even when we adopted a more restricted Toff-interval, and 
in that cases the procedure failed to isolate the best so- 
lution. In order to solve this situation, for PG 2131+066, 
PG 1707+427, and PG 0122+200 we restricted the search 
to within la of the spectroscopic value. 

We start by examining Fig. 1171 in which the run of 
is plotted versus Tcff for PG 1159-035. We include the 
results for the complete set of stellar masses considered 
in this work. The observed periods are the same twenty 
£ = l,m = consecutive periods considered in KB94. The 
vertical dotted line corresponds to the effective tempera- 
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Fig. 18. Same as Fig. [13 but for PG 2131+066. 



Table 3. Comparison between the observed periods of 
PG 2131+066 (taken from Kawaler et al. 1995) and the 
theoretical {£ — 1) periods of the best-fit model {SH — 1.83 
s, asu = 2.59 s). 







k 






[s] 


[s] 




[s] 


[%] 


341.45 


341.40 


14 


0.05 


0.014 




364.21 


15 






384.27 


382.23 


16 


2.04 


0.53 


403.93 


403.86 


17 


0.07 


0.017 


426.36 


425.07 


18 


1.29 


0.30 


450.28 


445.17 


19 


5.11 


1.13 


462.39 


466.26 


20 


-3.87 


-0.84 




486.74 


21 






507.91 


508.27 


22 


-0.36 


-0.07 



ture of PG 1159-035 as measured by means of spectro- 
scopic techniques {T^g°'^ = 140 000 K) . Note that there is 
a well-defined minima of x^, indicated with an arrow in the 
plot, corresponding to a stellar model with Af, — O.56M0 
and Toff ^ 137, 000 K. We adopt this model as our best-fit 
model. A comparison between the observed periods of PG 
1159-035 and the theoretical periods corresponding to the 
best- fit model is included in Table 2. The first column lists 
the observed periods, and the second and third columns 
correspond to the computed periods and the associated 
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radial order, respectively. The last two columns show the 
difference and the relative difference (in percent) between 
the observed and computed periods, respectively, defined 

as 5u, = np - and = 100 x (nf - nT)/np. 

Note that the period match is excellent. To quantitatively 
measure the quality of the fit we compute the mean of the 
period differences, 511 and the root-mean-square residual, 
asu- We obtain SU = 1.79 s and asn = 2.21 s for the 
fit to PG 1159-035. The quality of our fit is compara- 
ble to that achieved by KB94 {SU ~ 1.19 s), although 
the KB94 fit is notoriously better. In this connection, it 
is interesting to note that these authors adopt Mhc and 
the surface composition as free adjustable parameters for 
the period fit, in addition to the stellar mass and the ef- 
fective temperature. Therefore, they are able to make a 
finner tuning of the observed period spectrum than our 
own. The structural parameters of our best-fit model are 
listed in Table 6. Note the rather nice agreement between 
our predictions and the structural properties of PG 1159- 
035 as inferred by other standard techniques. For instance, 
the effective temperature, the luminosity and the surface 
gravity of our best-fit model are consistent at the la level 
with the spectroscopic estimations of Dreizler & Heber 
(1998). In particular, the quality function of our period fit 
naturally adopts an unique minimum value at an effective 
temperature very close to the spectroscopic value (this is 
not the case for the remainder GW Vir stars; see below). 
Note also that the total mass of our model lies at O.O2M0 
above the spectroscopic mass, although well within the 
allowed range (0.44 < M^/M© < 0.64). We stress again 
that spectroscopy determines only log 5 and Toff, and then 
the value of the stellar mass is assessed using evolution- 
ary tracks, which depend on the complete evolution of the 
progenitors. We note that our seismological mass is lower 
than the value obtained by KB94 by employing the same 
set of observed pulsation periods. 

At variance with PG 1159-035, the stars PG 2131-1-066, 
PG 1707+427 and PG 0122-1-200 exhibit very few pulsa- 
tion periods in their light curves. As we already antici- 
pated, this unfortunate fact renders any asteroseismic in- 
ference on these stars considerably more difficult. This is 
the reason for which most of asteroseismic studies carried 
so far rely on the measured period spacing only. In spite 
of this fact, we have attempted to apply the fitting 
method to these stars by adopting a more strict criteria, 
i.e., by considering as valid only solutions that guaranteed 
consistency at la level with the spectroscopic determina- 
tion of Teff. Fig. UHl shows the results for PG 2131+066. 
For this star the observational data consists of only the 
seven pulsation periods reported by Kawaler et al. (1995). 
Note that, at variance with the case of PG 1159-035, the 
function shows various minima compatible with T^g'^'^ 
in the range of effective temperatures considered. For in- 
stance, for = O.GOMq there is a remarkable minima 
in the immediate vicinity of the spectroscopic estimation, 
at Toff « 95 000 K. Other comparable minimum of 
are seen also for stellar masses of 0.61,0.63 and O.64M0, 
but at effective temperatures somewhat departed from 



Table 4. Comparison between the observed periods of 
PG 1707+427 (taken from Kawaler et al. 2004) and the 
theoretical {£ = 1) periods of the best-fit model {SU ~ 2.31 
s, asn = 2.70 s). 
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0.52 
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494.39 


489.62 


19 
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0.96 




514.23 
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536.41 


535.84 


21 


0.57 


0.11 




559.94 
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581.88 
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606.38 
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628.83 


25 








651.41 


26 






677.89 


677.48 


27 


0.41 


0.06 




697.36 


28 






726.02 


723.05 


29 


2.97 


0.41 


745.78 


747.18 


30 


-1.4 


-0.19 




769.01 


31 








792.73 


32 








818.54 
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841.27 


34 








862.57 


35 








890.30 


36 






909.05 


912.62 


37 


-3.57 


-0.39 



the spectroscopic value. These minima represent equiv- 
alent solutions from the point of view of the magnitude 
of x^, and thus, there are various possible optimal mod- 
els. We elect to adopt as the best-fit model to that having 
the effective temperature closer to the spectroscopic de- 
termination, that is, the model with Af» — 0.60Mq and 
Toff = 94 740 K (indicated with an arrow in Fig. HSJ). A 
comparison between the observed periods of PG 2131+066 
and the theoretical (^ = 1) periods of the best-fit model is 
given in Table 3. We note that the overall quality of the 
fit {SU — 1.83 s, asn = 2.59 s) is nearly so satisfactory as 
in the case of PG 1159-035 (Table 2). The characteristics 
of our best-fit model for PG 2131+066 are listed in Table 
6. Our results are in good agreement with spectroscopy, 
although our determination of the stellar mass is signifi- 
cantly higher — but even consistent at the la level. Note 
that the disagreement in the estimation of the stellar mass 
would be more pronounced if we were adopting anyone of 
the other acceptable solutions, for which > O.GOAf©. 

We have also applied our fitting procedure to the GW 
Vir stars PG 1707+427 and PG 0122+200. We have found 
that the behavior of in these cases is, not surprisingly, 
very similar to that shown by Fig.UHlfor PG 02131+066, 
because these stars also exhibit a very reduced number 
of pulsation periods. Therefore, to obtain optimal rep- 
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Table 6. Summary of the seismic inferences for PG 1159-035, PG 2131+066, PG 1707+427 and PG 0122+200 from 
period fits. In the interest of comparison, we include values of structural properties extracted from spectroscopic 
determinations and the values of the observed mean period spacing extracted from Table 6 of Kawaler et al. (2004). 



Quantity 


Model PG 1159-035 


Model PG 2131+066 


Model PG 1707+427 


Model PG 0122+200 


An [s] 

Teff [K] 

M^/Mq 
logs' [cm/s^] 
log(L,/Lo) 

R*/Re) 

Mae [10-=^ Mq] 
distance [pc] 


21.21 21.5 + 0.1 
137 620 140 000 + 6 000 
0.56 0.54 + 0.1 
7.31 7.0 + 0.3 
2.38 2.7 + 0.4 
0.0273 — 
4.94 — 
482 ± 44 800^40°'^^ 


20.86 21.6 + 0.4 
94 740 95 000 + 5 000 
0.60 0.55 + 0.1 
7.71 7.5 ± 0.3 
1.37 1.6 + 0.3 
0.0179 — 
5.29 — 


23.22 23.0 + 0.3 
87 585 85 000 + 4 500 
0.55 0.54 + 0.1 
7.59 7.5 ± 0.3 
1.31 1.4 + 0.3 
0.02 — 
4.85 — 


20.92 21.10 + 0.4 
81810 80 000 + 4 000 
0.64 0.54 + 0.1 
7.89 7.5 + 0.3 
0.96 1.3 + 0.3 
0.015 — 
5.65 — 



References: (a) Werner et al. (1991); (b) Reed et al. (2000) 



Table 5. Comparison between the observed periods of 
PG 0122+200 (taken from Vauclair et al. 2001)_and the 
theoretical {£ — 1) periods of the best-fit model {511 = 2.61 
s, cr^n = 3.55 s). 
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381.16 
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-1.56 


-0.41 


401.00 
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17 


-0.75 


-0.19 




424.05 


18 






449.43 


442.62 
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6.81 


1.51 


468.69 


463.65 


20 


5.04 


1.08 




485.45 


21 








506.63 


22 








527.22 
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546.99 


24 






570.00 


570.35 


25 


-0.35 


-0.06 




591.06 


26 






611.05 


611.24 


27 


-0.19 


-0.03 



rcscntative models wc employ the same procedure than 
in the case of PG 02131+066, i.e., we consider as valid 
the solutions that arc consistent at la with the spectro- 
scopically inferred effective temperature. In this way we 
discard other possible solutions at effective temperatures 
that do not match the spectroscopic evidence. The results 
are summarized in Tables 4 and 5, respectively. Note that 
for PG 1707+427 the period match is of a slightly lower 
quality than for PG_1159-035 and PG 2131+066, but stiU 
satisfactory, with 511 = 2.31 s and asn = 2.70 s. For PG 
0122+200, however, the fit is considerably worse, being 
511 = 2.61 s and asn = 3.55 s. The main source of dis- 
crepancy comes from the existence of periods at 449.43 
and 468.69 s. We mention that for PG 0122+200 we have 
also considered the set of observed periods reported by 
O'Brien et al. (1998), which do not includes the period at 
468.69 s. We find that the quality of the fit in that case do 
not improves significantly. The structural properties of the 
best-fit models for PG 1707+427 and PG 0122+200 are 
given in Table 6. For PG 1707+427 we found a general 
agreement between our inferences and the spectroscopic 



values, in particular concerning the stellar mass. In the 
case of PG 0122+200 we found a total mass and a surface 
gravity large in excess as compared with the spectroscopic 
evidence. This high mass value is in line with other seismic 
determinations and also with our own predictions based 
on the asymptotic period spacing and the average of the 
computed period spacings (see Tablel). 

In addition to the structural properties of the GW Vir 
stars undc;r study, we infer their "seismic" distance from 
the Earth. First, we compute the bolometric magnitude 
from the luminosity of the best-fit model, by means of 
Mboi = M0boi-2.51og(L*/L0), with Moboi = 4.75 (Allen 
1973). Next, we transform the bolometric magnitude into 
the absolute magnitude, Mv = Afboi — BC, where BC is 
the bolometric correction. Finally, wc compute the seis- 
mic distance according to \ogd [pc] = ^ (my — + 5). 
For PG 1159-035 (tov = 14.84) we adopt a bolometric 
correction of BC = —7.6 + 0.2 from KB94. Following 
Kawaler et al. (1995), we adopt BC = -6.0 ± 0.5 for PG 
2131+066 (mv = 16.6) and PG 1707+427 (rriv = 16.69). 
Unfortunately, there is not available any estimation of the 
bolometric correction for PG 0112+200, hindering any at- 
tempt to infer its seismic distance. Our seismic distances 
are shown in Table 6. In addition, we include the dis- 
tances estimated by means of other techniques. We found 
a good agreement (within Itr) between our distances and 
the other non-seismic estimations, although our values are 
characterized by errors much smaller. Our results are also 
consistent with the distances obtained seismologically by 
KB94 for PG 1159-035 {d = 440 + 40 pc) and by Reed et 
al. (2000) for PG 2131+066 {d = 668tg pc). 

5. Conclusions 

In this paper we have studied in detail some relevant as- 
pects of the adiabatic pulsations of GW Vir stars by using 
state-of-the-art evolutionary PG 1159 models recently pre- 
sented by Althaus et al. (2005). As far as we are aware, 
this is the first time that the adiabatic pulsation properties 
of fully evolutionary post born-again PG 1159 models like 
the ones employed in this work are assessed. We refer the 
interested reader to the paper by Gautschy et al. (2005) 
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for details about the non-adiabatic pulsation properties of 
these models. 

We first explored the basic nature of the pulsation 
modes by employing propagation diagrams and weight 
functions. In line with previous works, our results suggest 
that for PG 1159 models at high luminosity stages, the 
propagation diagrams are reminiscent of those of their pro- 
genitors — red giant stars — characterized by large values 
of the Brunt- Vaisala frequency in the central regions of the 
star. As a result, pulsation (7-modes — the only observed 
so far in GW Vir stars — are strongly confined to the 
highly condensed core, whereas p-modcs arc free to oscil- 
late in more external regions. In addition, we have found 
that there exist several modes exhibiting a mixed nature, 
behaving as p-modes in outer regions and like g-modes in 
the deeper zones within the star. As the effective temper- 
ature increases these modes undergoes several episodes of 
avoided crossing. We found that pulsation periods gener- 
ally decrease with time, an effect attributable to the rapid 
contraction experienced by the star in their incursion to 
the blue in the H-R diagram. Note that a decreasing period 
implies negative temporal rate of period change; accord- 
ing to our results, this is the case for all the modes during 
the contraction stage. Once the models have passed their 
maximum effective temperature and have settled onto the 
WD cooling track, the Brunt- Vaisala frequency acquires a 
more familiar shape, typical of the WD pulsators. In this 
phase all the pulsation periods increase with decreasing 
effective temperature, refiecting a lowering in the magni- 
tude of the Brunt- Vaisala frequency at the core regions. 
At this stages, (7-modes become envelope modes, and, as 
indicated by the weight functions, the outer regions of the 
star are the more relevant ones in establishing the periods. 
By the contrary, pulsation p-modes are almost confined to 
the degenerate core. We note that, since all the periods in- 
crease during this phase, the rate of period change for any 
pulsation mode must be positive. In particular, the mea- 
sured rate of period change of the period at ~ 516 s in 
PG 1159-035 is positive (Costa et al. 1999), in agreement 
with our findings. We defer to a future work a complete 
discussion of the rate of period change of our PG 1159 
models. 

We next have focused our attention on the evolution- 
ary stages in which the star has already entered its WD 
cooling track. For these phases we have been able to ob- 
tain additional evolutionary sequences with several values 
of the stellar mass, allowing us to study the effects of M» 
and the effective temperature on the pulsation properties 
of our models. In particular, we have examined the asymp- 
totic behavior of the (7-modc; pulsations. In agreement with 
previous works, we found that the asymptotic period spac- 
ing increases with a decrease in the stellar mass and with 
an increment in the luminosity, although the dependence 
on the stellar mass is stronger. This fact renders the value 
of the period spacing observed in a given star a powerful 
indicator of the total mass. 

The study of the mode- trapping properties of our mod- 
els has been other relevant aim of the present work. In 



this connection, we have demonstrated that the mode- 
trapping features characterizing our PG1159 models are 
produced mostly by the shape of the 0/C chemical pro- 
file at the core, at least for the range of periods observed 
in GW Vir stars. On their hand, the outer chemical in- 
terface of 0/C/He produces negligible mode-trapping ef- 
fects, except for the regime of short periods in the pulsa- 
tion spectrum. This conclusion is at odds with previous 
results reported by KB94, who realized that the mode- 
trapping properties of their PG 1159 models were fixed 
essentially by the outer 0/G/He transition region. We 
have found that the origin of this discrepancy rests on 
the differences in the input physics employed in build up 
the background stellar models. Specifically, the main dif- 
ference is the presence of a much less pronounced chemical 
transition in the outer parts of the C/0 core of the KB94 
models, as compared with the rather abrupt chemical gra- 
dients at Mr « 0.4-0.6M* in our PG 1159 models. Other 
ingredient that contributes to the disagreement between 
our results and those of KB94 is the thickness of the he- 
lium envelope. Indeed, since our models are characterized 
by thick helium envelopes (Mne ~ O.OO52M0), the mode 
trapping effects caused by the 0/C/He transition region 
are much more weaker than in the models of KB94. 

We note that the structure of the core chemical profile 
of our PG 1159 models is the relic of convective over- 
shoot episodes during the central helium burning phase. 
The sensitivity of the mode trapping effects to the de- 
tails of the core chemical structure rises the possibility 
of employing pulsating PG 1159 and WD stars to con- 
strain the efficiency of extra mixing episodes in the core 
of their progenitors (overshooting and/or semiconvection; 
see Straniero et al. 2003). This appealing issue has been 
recently explored by Corsico & Althaus (2005). 

To gain additional insight into the nature of mode 
trapping in our models, we have examined the kinetic en- 
ergy and weight functions of individual pulsation modes. 
We have found that modes showing primary maxima in 
the kinetic energy distribution are associated with min- 
ima in the period-spacing diagrams. By consulting the 
cigcnfunctions and weight functions of these modes we 
found that they have relatively high amplitude at the 
high-density environment of the stellar core, and very low 
amplitudes in the rest of the star. The existence of these 
"core-trapped" modes is encountered also in the context 
of ZZ Ceti pulsations, as reported by Althaus et al. (2003) 
and Corsico et al. (2005). 

Finally, we have made some preliminar seismic in- 
ferences about the internal structure and basic parame- 
ters of the GW Vir stars PG 1159-035, PG 2131-^066, 
PG 1707-F427, and PG 0122-^200. To this end we have 
adopted three different approaches. First, we estimate the 
stellar masses by comparing the asymptotic period spac- 
ing with the observed mean period spacing, assuming that 
the effective temperature of the target star is that pre- 
dicted by spectroscopy. In the second method, we repeated 
the above procedure but this time using the average of 
the computed period spacings and comparing it with the 
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observed mean period spacing. The third approach is a 
fitting method in which we search for an optimal stellar 
model that best reproduces the observed periods. To this 
end we employ a quality function that measures the dis- 
tance between the observed and the computed adiabatic 
pulsation periods within a grid of models with different 
effective temperatures and stellar masses. 

The stellar masses we obtained from these methods 
are given in Table 1. Note that, except for PG 0122+200, 
our values are consistent with the spectroscopic inferences, 
particularly when the average of the computed period 
spacings are used, although our masses are slightly larger. 
For PG 0122+200 we obtain, in agreement with the seis- 
mic determination of O'Brien (2000), a large value of the 
stellar mass. Thus, according the pulsation theory this star 
seems to be the most massive among the high-gravity GW 
Vir pulsators. Note also that the stellar masses as deter- 
mined by employing the average of the computed period 
spacings are appreciably lower than the values inferred by 
using the asymptotic period spacing. We also found that 
our M, values are smaller than that reported by other 
pulsation studies. In this connection, we note that, with 
the exception of the work of KB94 for PG 1159-035, all 
these works compare the observed mean period spacing 
with the asymptotic predictions to infer the stellar mass. 
Thus, not surprisingly, their values are very close to our 
results obtained from the asymptotic period spacing, but 
somewhat different from our values deduced by employing 
the average of the computed period spacings. 

The results from our period-fitting procedure are 
shown in Table 6. Taking full advantage of the numer- 
ous pulsation periods observed in PG 1159-035, we have 
unambiguously obtained a best-fit model with an effective 
temperature very close to that predicted by spectroscopy. 
Our asteroseismic mass of O.56M0 is consistent with the 
spectroscopic calibration (Af* = O.54M0) and with the 
preferred value of Gautschy et al. (2005) (M* = O.55M0), 
but considerably lower than the value quoted by KB94 
{M^, = O.59M0). For the remainder stars, for which we 
have very few observed periods available, we have also 
obtained representative models by employing our period- 
fitting procedure, but, at variance with the case of PG 
1159-035, for these stars we have been forced to consider 
as acceptable solutions only models with an effective tem- 
perature laying within la of the spectroscopic value. This 
stringent criteria was necessary in view of the numerous 
and almost equivalent minima exhibited by the quality 
function for these stars. Note that the period matching 
are considerably poorer (Tables 3, 4 and 5) as compared 
with the case of PG 1159-035 (Table 2). We also have es- 
timated the "seismic" distances to the stars by employing 
the luminosity of the best-fit models. Our estimations (see 
Table 6) are consistent with other determinations. 

Before closing the paper we stress that, in our view, 
the pulsational analysis presented here constitutes a sub- 
stantial improvement as compared with previous studies. 
However, we believe that PG1159 evolutionary models 
with different stellar masses based on a complete descrip- 



tion of the physical processes occurred in all the evolution- 
ary stages of progenitor stars would be needed to reinforce 
some of our results. The development of such models is cer- 
tainly critical for the assessment of the internal chemical 
composition characterizing GW Vir stars, an issue that is 
key as far as precision asteroseismology is concerned. 

Finally, the evolutionary tracks employed in 
this investigation, as well as tabulations of Hq 
in terms of L* and T^ff for the complete set of 
stellar masses are freely available at our URL: 
|http : / Twww . f caglp . unlp . edu . ar/evolgroup/ 
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Appendix A: Pulsation equations and adiabatic 
quantities 

Our numerical pulsation code solves the fourth-order set 
of equations governing linear, nonradial, adiabatic stellar 
pulsations in the formulation given in Unno et al. (1989): 



dyi 
dr 

dy2 
dr 

dr 

dyi 
dr 



2/2 + 2/3, (A.l) 



- A*) + (A* - [/ + 1) y2 - A* ys, (A.2) 
(1 - U) j/3 + y4, (A.3) 



UA* yi+UVg y2 + W + 1) - UVg] ys-U 2/4. (A.4) 
The boundary conditions are, at the stellar centre (x — 0): 
y, CiLo^-£y2 = 0, £ ys - ^ 0, (A.5) 

and at the stellar surface {x = 1): 

2/1 - 2/2 + ys = 0, {£ + 1) + 2/4 = 0, (A.6) 

being the normalization yi = 1 at a: = 1 (x = r/R^). 
The dimensionlcss Dzicmbowski's variables (eigenvalue 
and eigenfunctions) are defined as 



and 



2/1 



2/2 



2/3 = — 
gr 



1 

2/4 = --r- 

9 dr 



(A.7) 



(A.8 



(A.9) 



Here, £,r is the radial Lagrangian displacement, and p' , 
<&' are the Eulerian perturbation of the pressure and the 
gravitational potential, respectively. Pertinent dimension- 
less coefficients of the pulsation equations are: 



V 



u = 



Airpr^ 

Mr 

A* = 



9 



(A.IO) 



(A.ll) 



where c is the adiabatic local sound speed and N the 
Brunt- Vaisala frequency. The remainder symbols are self- 
explanatory. Once the eigenvalue and eigenfunctions are 
computed, the code proceeds to evaluate a number of im- 
portant pulsation quantities, such as the pulsation period, 



Uk - 2n/ak, 

the oscillation kinetic energy, Kk, 



(A.12) 
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Kk = -{GM.Ri)u;i 



X / X p 



2 2 I 2_^(^_j_J_) 2 

X yi -\- X 1^ 2N2 yi 



the weight function, Wk, 



dx, 



(A.13) 



A*y\ + V, (z/2 - 2/3)' - ^{^ {I +1)2/3 + 2/4}' 



the variational period, 11^, 
Jo 



1/2 r „1 



-1/2 



(A.14) 



(A.15) 



and finally, the first-order rotation splitting coefficients, 
Cfe, 



2a; 2/12/2 



2/2 



dx, (A.16) 



The rotation splitting of the eigenfrequencies (assuming 
slow, rigid rotation) can be assesed by means of 



CTfe.TO = cTfe + m(l - Cfe)0 



(A.17) 



where O is the angular speed of rotation and m the az- 
imutal quantum number. 



